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Abstract 

The work of formation of a critical nucleus is sometimes written as W = 
nAfi + 7 A The first term W vo i = nAfi is called the volume term and the 
second term jA the surface term with 7 being the interfacial tension and A 
the area of the nucleus. Nishioka and Kusaka [J. Chem. Phys. 96 (1992) 
5370] derived W vo i = nAfi with n = Vp/vp and Afi = fip(T ', p a ) — fi a (T ', p a ) by 
rewriting W vo i = —{jpp — p a )V/3 by integrating the isothermal Gibbs-Duhem 
relation for an incompressible (3 phase, where a and (3 represent the parent 
and nucleating phases, Vg is the volume of the nucleus, v@, which is constant, 
the molecular volume of the (3 phase, fi, T, and p denote the chemical po- 
tential, the temperature, and the pressure, respectively. We note here that 
Afi = fipiT^pa) — fi a {T,Pa) is, in general, not a directly measurable quantity. 
In this paper, we have rewritten W vo i = —(pp — p a )V/3 in terms of \i re — /i eg , 
where /i re and fi eq are the chemical potential of the reservoir (equaling that 
of the real system, common to the a and phases) and that at equilibrium. 
Here, the quantity \i re — fi eq is the directly measurable supersaturation. The 
obtained form is similar to but slightly different from W vo i = nAfi. 
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1. Introduction 



To calculate the reversible work of formation of a critical nucleus is one 
of the purposes of the theory of nucleation, because one can predict the 
steady state nucleation rate J s = J exp(—W*/k B T) through the work of 
formation of the critical nucleus W, where k B T is the temperature multiplied 
by Boltzmann's constant. Here, W* = W(R*) is the height of the nucleation 
barrier with R* being the radius of the critical nucleus. We can refer a theory 
not including molecular level quantities to as a classical nucleation theory. 
We often encounter the following formula [Eq. (CQ)] or equivalent one in the 
classical nucleation theory: 

W = nAfi + ^A, (1) 

with 7 begin the interfacial tension, A = AirR 2 the area of the interface 
(rigorously speaking, R is the radius of the surface of tension) in textbooks 
such as (H0jyj0jSjla| as well as research papers such as (7l.l8l.l9l.E0. 11, 12. 13 



Nishioka and Kusaka 14] called this formula "a commonly used formula." 
In this formula, one regards A/i the chemical potential difference between 
the parent phase (the a phase) and the nucleating phase (the j5 phase), i.e., 
Afj, = up — fi a (< 0) [in the field of crystal growth, one sometimes defines the 
supersaturation as A/x = /x a — fAp(> 0) and the negative sign arises as — nA/i]. 
We note here that the term "supersaturation" is sometimes used to express 
the such thermodynamic driving force, apart form the literal meaning. In 
this paper, we will adopt this terminology. One can understand n as the 
numbers of molecules undergone the phase transition from the a phase to 
the (3 phase. Apart form its implication, meaning of A/i in this expression 
is immediately ambiguous; for the critical nucleus, the chemical potential 
of the (3 phase is equal to that of the a phase. In this paper, we will give 
correct formulas for W in terms of the chemical potential difference relative 
to equilibrium one. 

is 



The exactly correct form for W given by Gibbs 15 



W = -( Pf} -p a )Vp + yA, (2) 

where p denotes the pressure and Vp = 47i\R 3 /3 is the volume of the nucleus. 
Rigorously speaking, pp is the pressure of the hypothetical cluster defined 
such as possessing the bulk property and filling the inside of the surface of 
tension, and thus Vr the volume inside the surface of tension. A transparent 
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explanation for the volume term W vo \ = — (pp — p a )Vp can be given through 
a grand potential formalism — the grand potential is defined as Q = —pV; 
the grand potential formalism for interfaces was given in, for example, a 
textbook by Landau and Lifshitz [l6[ and exactness of W vo i — ~{Pp ~ Pa)Vp 



is shown by a textbook by Vehkamaki 17] (see also 18|). The author wish 
to introduce a heuristic paper [19| for readers' convenience. That is, one 
can readily understand the form of W vo \ on the basis of the fact that the 
reversible work of formation of the critical nucleus is the grand potential 
difference. In the following sense, Eq. (j2J) is entirely exact. We divide the 
process of nucleus formation into two. One is the formation of a hypothetical 
cluster of radius R, within which the bulk (3 phase fulfills. The other is the 
formation of interfacial structure on the mathematical boundary of radius R. 
The former is given by the first term in Eq. (j2J) and the latter is expressed 
by the second term. The problem of classical theories is that the interfacial 
tension 7 is treated as constant (in particular, in a capillary approximation 
7 for the flat interface is employed), whereas the curvature dependence is 



unnecessarily negligible [20|, |21|, 122, [23|, [24], l25|, |26|, |2 



In Sec. [2] we will give the derivation of W vo i = n[fip(T,p a ) — fi a {T,p a )] 

Here, we will give a brief review on this 
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due to Nishioka and Kusaka 
form. To express the pressures in terms of the chemical potentials using 
the Gibbs-Duhem relation was seen in Oxtoby and Kashchiev's paper on the 
nucleation theorem 
Laaksonen et al. 
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In their paper, however, W vo i was not calculated. 
291 ] pointed out that W vo i could be calculated by extending 
the Oxtoby and Kashchiev's line. This is completely the same as Nishioka 
and Kusaka' s 



and Reiss 3^ 
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Also, the same procedure was followed by Debenedetti 
In order to avoid defocusing of the point of the present paper, 
we will not give a general review on the nucleation theorem. 

Differentiating Eq. (J2J) with respect to R and equating zero and then 
solving for R, we have formula for the size of the critical nucleus as 



R* 



2 7 



Pp-Pc 



(3) 



By substituting Eq. ([3]) for R in Eq. (j2~J) the hight of the nucleation barrier 
is obtained as 



W* = -(pp - Pa )V$ = ->yA* 



167T7 3 



KPP -Pa) 2 ' 



(4) 



.3 



with Vg = 47r(i?*) 3 /3 and A* = 4n(R*) 2 . Corresponding to Eq. ([I]), one has 



2„,3 



W* 



167TV 7 

3(A^ 



(5) 



instead of the last expression of Eq. (J3J). Here, t> represent "the molecular 
volume"; that is, n (or n*) is expressed as Vp/v (or Vg/v). We note here that 
the definition of v is immediately ambiguous; whose phase is not specified, 
or v may be common to the a and /3 phases. Also, as mentioned above, the 
definition of A/i immediately unclear. For a rarefied a phase, one uses very 
often 

A/i = ksT In S, (6) 

where S is the supersaturation ratio p re /Peq with p re being the pressure of 
the reservoir, which equals p a , and p eq the equilibrium pressure. In other 
words, A/i is defined as 



A/i = /i a (T,p re ) - fi a (T,p 



eq)i 



(7) 



where T is the temperature, which is assumed to be uniform throughout the 
system. Hereafter, T will be omitted for brevity. As shall be shown below, 
this definition is, however, incorrect. Nevertheless, expression W* in terms 
of this A/i is strongly desired because this quantity is directly measurable. 
Indeed, Eq. ([7]) is widely used such as in [3l|, |32|, |33|, |34|, |35|, |36|, [37 . 



2. Derivation of commonly used formula 



Nishioka and Kusaka 14] found out that in case that the (3 phase is 
incompressible, the volume term W vo i = —{pp — p a )Vp can be rewritten in 
the form n(p,p — fi a ). That is, they integrated 

V ' - (8) 



dp 

which is nothing other than Gibbs-Duhem relation for the isothermal case, 
for the P phase. Unfortunately, they concluded incorrectly that the form 
([1]) was valid only for the case of the incompressible /3 phase such as the 
nucleation of an incompressible liquid phase in a vapor phase. For example, 
this condition is, in a mathematical form, valid for a bubble nucleation in 
an incompressible liquid phase. In this paper, however, we limit ourselves to 
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the incompressible (3 phase to avoid the confusion in argument. Note that 
the form of Eq. (pQ) cannot be derived in general — we can derive in some 
approximations. 



For the case of the incompressible (3 phase, W vo i is given by [14 



Wvol = —\pLp{p a ) - Haip a )\. (9) 
Vfi 

To reach to this expression we have used the fact that the chemical potential 
of the nucleation phase Hpipp) is equal to that of the parent phase /j, a (p a ), 
i.e., 

PpM = Va(Pa), ( 10 ) 

(that is, the chemical potential is uniform throughout the system) as indi- 
cated by a horizontal dashed lines in Fig. [TJ The vertical solid lines in Fig. [1] 
depicts Hp{p a ) — fJ>a(p a ) [Eq. OH])]. Figured] (a) is for a normal case such as 
a liquid droplet in a vapor phase. On the other hand, Fig. [1] (b) is for an 
abnormal case such as formation a nucleus of less denser, incompressible /3 
phase in a denser a phase. In both cases, np(p a ) —^a{Va) does not coincide to 
the supersaturation, which is the chemical potential difference relative to the 
chemical potential at the a-/3 phase equilibrium (sometimes with a negative 
sign as note in Sec. [T]). 



3. Volume term in general 

In these ways, one knows two cases where the form of eq. (pQ) is valid. Let 
us develop a general consideration; we consider the case that v in eq. (JSJ) is a 
function of p according to the mean value theorem. Following Nishioka and 
Kusaka [l4| we integrate the Gibbs-Duhem relation of the form of eq. (151) 
from pi to p2- 

KP2) ~MPi) = v (P)(P2-Pi), (11) 

where p is a certain value lying in the interval (p2,Pi)- Let us define p eq as the 
equilibrium pressure of the a-(3 phase equilibrium (the saturation pressure 
of the a phase with respect to the j3 phase), which is the solution to 

t*p(p) = Va(p)- (12) 

For latter convenience, let us define an abbreviation fi eq for Hp(peq) = ^a(p eq )- 
Also for latter convenience, we define /i re = Hpijip) = A*a(Pa)> which is the 
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chemical potential of the reservoir in /iVT ensemble. Applying eq. ( II ip for 
the (5 and a phase, we have /ip(pp)- Hp{peq) = vp{p eq + 6p(pp-p eq ))(pp-p eq ) 
and Ha(p a ) - Ha(p eq ) = v a (p eq + 9 a [p a - p eq ) ) (p a - p eq ) with Op and 9 a being 
certain values laying in an interval (0,1), i.e., 

Hre-Veq = Vp(fp) (Pj8 ~ Peq) , (13) 
flre-^eq = V a (f a ) (j) a - p eq ) , (14) 

where and t> Q are molecular volumes of respective phase, p~p and p~ a certain 
values respectively lying in intervals {p eq ,Pp) and (p e q,p a ) for cases of pp > 
Pa > Peq as shown in Fig. [2] (a) [for cases of p a < pp < p eq as shown in 
Fig. 121(b) the intervals are replaced with (pp, p eq ) and (p a , p eq ), respectively]. 
By dividing Eq. ( Tl3|) by vp{p~p) and Eq. (Tl4j) by v a {p a ) and subtracting the 
latter from the former, we eliminate p eq to have 
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(/Ire ~ Heq) = P/3 ~ Pa- (15) 



Vp(pp) v a (p c 

We can rewrite the first term in Eq. (|2l) by substituting pp — p a by Eq. (US]). 



(/' re f^eq )Vp 



= -{rip - n tt )(/i re - fi eq ), (16) 

where np = V^/u/3(p~/3) and n a = Vp/v a (p a ). This form is slightly different 
from the first term in Eq. ([1]). In a case that the a phase is a rarefied 
vapor phase, the quantity n a tends to vanish and the form of the first term 
in eq. (jTj) is obtained. This case can be included in the case of Fig. [TJ (a) 
(nucleation of an incompressible liquid droplet in an infinitely rarefied vapor). 
In this case the curve of \x a tends to a vertical line. In this limiting case the 
horizontal location of the intersection of horizontal dashed line with the fi a 
curve locates on the intersection of two \i curves; that is, fip(p a ) — Ha{p a ) in 
eq. ([H]) tends to coincide to /i re — fi eq . Here, we note that — — fi eq ) is the 
true supersaturation for the case of Pp > p a > p eq [Fig. [21 (a)] and the true 
undersaturation for the case of pp < p a < p eq [Fig. El (b)]. In this way, we 
have successfully express the work term in terms of the supersaturation and 
revealed under what case the commonly used formula holds. 

We have obtained a form of Eq. ([T]) with A/x = /i re — fi eq . The result 
of Eq. (fT5"|) includes an issue concerning physicochemical problems, although 
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it is of a mathematically beautiful form; in definitions rip = Vp/vp(p~p) and 
n a = Vp / v a (p a ) , the denominators vp(p~p) and v a (p a ) are not fixed constant 
values. We will give a prompt solution here. Assuming a smallness of \ire~Peq 
[rigorously speaking the smallness should be described in terms of a dimen- 
sionless quantity — that is, the present statement reads |(/z re — p eg )/ p eq \ <C 1 
or \fi re — fi eq \/k B T <C 1], let us make an second-order expansion instead of 
Eqs. ( I13p and (TH|) . Instead of Eq. OH]), it is more convenient to start with 

=A (17) 



dfJLy T 

with p = 1/v being the number density. The second-order expansions are 

1 



P/3 ~ Peg = pp(p e q)(Pre ~ Peq) + \p\kV 'eq) K fi(p 'eq) '(P re ~ Peqf 

3^ 



+0((/i re -/i e9 ) 3 ), (18) 

Pa - Peq = Pa{Peq){Pre ~ Peq) + ^ Pa (Peg) K a {Peq) {.Pre ~ Peqf 

+0((p re -p eq ) 3 ), (19) 

where k denotes the isothermal compressibility. Subtraction Eq. ffl9l) from 
Eq. ffl8|) . we have 



Pp-Pa = \pp{Peq) ~ Pa{Peq)]{Pre ~ Peq) 

+ ^[pl{Peq)Kp{Peq) ~ P a {Peq) ^ (Peq)] {Pre ~ Peqf 

+0{{p re - p eq f). (20) 

Neglecting the second and higher order terms in Eq. ( 1201) and inserting in 
the first term of Eq. (T5]), we have the form of Eq. ( TloT) with fixed rip and n a 
as an approximation. Taking into account the second term, we can improve 
the approximation. 

Coincidence between the present result and the commonly used formula 
with Eq. (J7j) can be understood as follows. For nucleations of incompressible 
f3 phase in a rarefied gas pp(p a ) — PaiPa) can be rewritten as 

PpiPa) ~ PaiPa) = [PpiPa) ~ Peq] ~ [Pa{Pa) ~ Peq] 

= MPa) - Pp(Peq)} - [Pa(Pa) - Pa(Peq)} 

= vp(p a - p^) - k B T In— . (21) 

Peq 
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The last term can be expanded in p a — p eq as 



/c^Tln — 

Peg 



ksT In 



Pc 



1 + 

-Peg 



Pa-P 



<<! 



P 



eg 



k B T 

Peg 

Va{Peg){p a P 



eg ) 



(22) 



To reach to the last line, the equation of state for the ideal gas has been used. 
Because the a phase is a rarefied gas, the inequality v a 3> Vp holds and then 
we find that the first term in Eq. (12"T]) can be neglected. 



4. Discussion 



In 1984, Wilemski 38] divided the number of molecule of species i in- 
cluded in a nucleus for the binary system as 



n b + nf, 



(23) 



with the superscripts b and s denoting bulk and surface and wrote down the 
condition of the critical nucleus as 



= (A/i + 7(aA/n i ) n J(d7ij)n i 
+n b 1 d/j,[ + n b 2 d/J 2 

+n{dii l 1 + n s 2 dfi l 2 + Ad>y (const. T,P). (24) 

In 1999, Laaksonen et al. [29| revealed the work of formation of the nucleus 
underlaying Eq. ( |24"|) as 

AG = ^(M^) - Hvi(Pv))9i + A-y, (25) 

i 

with gi = Tin — Ti v i + n si , where the subscripts 1 and v denote the liquid 
and vapor phases (in Wilemski's paper, the starting equation is AG = 
TT-i A/xi + ri2A[i 2 + Ay). The so-called surface excess number of molecules, 
n s , is identical to the superficial number of molecules in the Gibbs interfacial 
thermodynamics [l5| and we find that n s is proportional to A. 

As mentioned in Sec.[H in the Gibbs interfacial thermodynamics the work 
of nucleus formation is divided into the formation of the hypothetical cluster 
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and that of the interfacial structure. A quantity proportional to a superficial 
quantity is categorized as the latter. In this respect, in Eq. (|2"5|) the volume 
term is regarded as 

W vo i = y~](jJ.u{Pv) - fJ*vi(Pv)){nu - n vi ) (26) 

i 

While the factors regarded as A/i are different with each other in Eqs. fflBT) 
and (126]) . the coefficients to "A/i" coincide with each other. The difference 
is that in Eq. ff25l) the quantity fJ,p(p a ) — fJ>a(jPa) appears, not the quantity 

5. Concluding remarks 

We have successfully rewritten the volume term of the work of formation 
of a critical nucleus in terms of the supersaturation. The result is similar 
to the form of W vo i = nA/i but slightly different; n in this form has been 
replaced with np — n a . The form W vo i = nA/i with A/i being the supersatu- 
ration (the chemical potential difference relative to equilibrium) is recovered 
in a limiting case that the parent phase is a rarefied gas. This is a find- 
ing that requires a concerning note in textbooks — as mentioned in Sec. (U 
some textbooks lead readers to understanding that W vo i = nA/i is exact, and 
in some literatures this form is valid only for an incompressible nucleating 
phase. 

We wish to postpone comparisons with experimental studies after formu- 
lation in a form of the nucleation theorem. That is, W* should be plotted 
against experimentally determined A/i's to evaluate the deferential coeffi- 
cient. In relation with experiments, a crucial comment arises: the mathe- 
matical formula itself (without expansion made in the latter part of Sec. [3]) 
makes a sense — if one forcibly write the work of formation of a critical 
nucleus in the form like Eq. (JTj), then uncertainties necessarily accompany in 
interpreting the experimental results. 
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Figure 1: fi-p relation for cases of incompressible f3 phases; (a) a normal case that the (3 
phase is denser than the a phase and (b) an abnormal case that the a phase is denser than 
the /3 phase. The horizontal dashed lines indicate the chemical potential, which is common 
to the nucleus and the parent phase. The vertical solid lines depict np{p a ) — p. a (p a ) 
[eq. ©]. 




Figure 2: Location of /j eq , pp, p a , and p eq is illustrated (note that the reserver pressure p re 
coincides to p a ); (a) a case of condensation, i.e., the j3 phase is denser than the a phase, 
and (b) a case of bubble nucleation, i.e., the a phase is denser than the (3 phase. In case 
(a), p/3 > p a > Peg holds. On the other hand, in case (b), p a < Pp < Peg holds. 
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